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The occupation of more than one single-particle state and hence the emergence of fragmentation 
is a many-body phenomenon universal to systems of spatially confined interacting bosons. In the 
present study, we investigate the effect of the range of the interparticle interactions on the fragmenta¬ 
tion degree of one- and two-dimensional systems. We solve the full many-body Schrodinger equation 
of the system using the recursive implementation of the nrulticonfigurational time-dependent Hartree 
for bosons method, R-MCTDHB. The dependence of the degree of fragmentation on dimensionality, 
particle number, areal or line density and interaction strength is assessed. It is found that for con¬ 
tact interactions, the fragmentation is essentially density independent in two dimensions. However, 
fragmentation increasingly depends on density the more long-ranged the interactions become. The 
degree of fragmentation is increasing, keeping the particle number N fixed, when the density is 
decreasing as expected in one spatial dimension. We demonstrate that this remains, nontrivially, 
true also for long-range interactions in two spatial dimensions. We, finally, find that within our fully 
self-consistent approach, the fragmentation degree, to a good approximation, decreases universally 
as IW 1 / 2 when only N is varied. 

PACS numbers: 03.75.Gg, 03.75.Nt 


I. INTRODUCTION 

The many-body behavior of systems of interacting 
bosons can in many cases be essentially characterized 
by a single “diluteness parameter” expressing the im¬ 
portance of interactions |T]. For most experimentally 
produced samples of ultracold bosonic gases, describing 
them by such a (small) diluteness parameter essentially 
yields the correct physics, because the mean interparti¬ 
cle distance is much smaller than the range of the in¬ 
terparticle interaction potential. Archetypical theoreti¬ 
cal methods employing the diluteness of the gas are the 
Gross-Pitaevskii mean-field theory [T] and, with (quan¬ 
tum) fluctuations incorporated to lowest order, Bogoli- 
ubov theory [2]. These approaches assume that a single 
macroscopically occupied orbital forms the many-body 
ground state, and hence that a single creation operator 
in the field operator expansion is sufficient. 

In particular in one and in three spatial dimensions 
(ID and 3D), diluteness criteria are well established for 
integrable as well as nonintegrable interparticle interac¬ 
tions with a potential V[ n t(r) oc l/r n , [3]. Integrabil- 
ity here means the (infrared) convergence of the integral 
I Vi n t(r)d D r, implying n > D, where D is the spatial 
dimension. In the following, we use the terminology of 
interactions being short- and long-range interchangeably 
with them being integrable and nonintegrable, respec¬ 
tively. It is a widely accepted fact that for short-range 
interactions, the gas is more dilute and theories such as 
Bogoliubov mean-fields are applicable for low densities in 
3D and, conversely, for large densities in ID. For long- 
range interactions, that is in particular for Coulomb in¬ 
teractions, the gas cannot be considered as dilute when 


the density decreases in 3D. A consequence of this fact 
is, for example, the formation of a Wigner crystal [3|. 
The 3D case can however for long-range interactions like 
Coulomb be described by Bogoliubov theory for suffi¬ 
ciently high densities [5]. On the other hand, in one di¬ 
mension, high densities, in combination with Coulomb 
interactions, that is n = 1, lead to a Tonks-Girardeau 
gas 0- For strongly repulsive contact interactions in ID 
the density profile can either be captured by a Thomas 
Fermi or a Tonks-Girardeau local equation of state for 
high or low densities, respectively [7]. 

Two spatial dimensions (2D), which are the main fo¬ 
cus in what follows, represent a crossover case between 
the two limiting behaviors of ID and 3D with respect 
to the effective diluteness of the gas, and therefore con¬ 
stitute a sensitive probe for many-body effects. In 2D, 
correlations can potentially play a significant role for 
both, integrable and nonintegrable interactions. For in¬ 
tegrable interactions like contact pseudopotentials it is 
known that the diluteness parameter in 2D is essentially 
independent on the areal density p , involving a double 
logarithm In In [pa 2 ], where a s is the s-wave scattering 
length [S]. In the infinitely extended homogeneous gas, 
two dimensions represent the marginal case for the ex¬ 
istence of condensates EM; at any finite temperature 
condensation is absent, while at zero temperature con¬ 
densation is possible albeit for extremely small densities 
Il3l . Two-dimensional Bose gases have been comprehen¬ 
sively studied for short-range interactions, cf. the re¬ 
view M- However, comparatively little is known for 
long-range interactions, apart from the case where gauge 
fields of constant large flux are externally applied to two- 
dimensional electron gases, where the theoretical founda- 
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tion of the fractional quantum Hall effect in the form of 
the eponymous Laughlin wavefunction m has spurred 
both tremendous experimental and theoretical activity 
m- We also note in this regard that the many-body 
physics of dipolar bosons in 2D, cf., e.g., El does not 
qualify with the presently employed terminology as be¬ 
ing genuine “long-range,” due to the integrability of the 
interaction potential of (polarized) dipoles in two spatial 
dimensions. 

The question how the diluteness of a bosonic many- 
body system depends on the power law tail of the inter¬ 
action potential has thus not yet been addressed gen¬ 
erally, in particular for the two-dimensional case. In 
the present work, we aim to fill this gap by employing 
the recursive implementation of the multiconfigurational 
time-dependent Hartree for bosons (R-MCTDHB) !18j . 
to solve the many-body problem with both integrable and 
nonintegrable interaction potentials of varying power law 
n. 

The phenomenon of fragmentation, cf., e.g., [TTi! - f2TI] ■ 
in which more than one field operator mode becomes 
macroscopically occupied, so that Bogoliubov mean-field 
theory necessarily breaks down, has previously in par¬ 
ticular been extensively studied in one dimension. For 
instance, a stepwise increase of the number of fragments 
takes place in an asymmetric double well [21]. Because 
of the strongly increased demand in computational re¬ 
sources the two-dimensional and three-dimensional cases 
have only recently attracted attention due to the avail¬ 
ability of novel numerical methods, see for instance El- 
Fragmentation has been found in the case of an effective 
potential barrier [23] induced by the long-range interac¬ 
tions or the dynamics above a condensed ground state 

m- 

We focus below on studying fragmentation as an “in¬ 
trinsic” property of the many-body system in a single 
(e.g. hard wall or harmonic) trap. In this case, fragmen¬ 
tation is termed intrinsic, because it is not assisted by 
the one-body potential in the Hamiltonian. This is in 
contrast to the cases where fragmentation is “extrinsic”, 
i.e., assisted by the one-body potential as for instance 
in double wells or optical lattices. We refer to intrinsic 
fragmentation as a genuine many-body effect, because 
it originates solely in the interparticle interactions and 
thus does not essentially depend on the way in which the 
bosons are confined. 

Exactly solvable cases for benchmarking the numer¬ 
ics are very scarce and significantly more specialized 
in spatial dimensions higher than one. For the one¬ 
dimensional case, exactly solvable cases are available such 
as the Tonks-Girardeau gas and the Lieb-Liniger and 
Calogero-Sutherland models, as reviewed in 1223 ■ Analyt¬ 
ical ground state solutions in two dimensions have been 
obtained for instance for the bosonic variant of Laugh- 
lin’s wavefunction for particles in an (effective) magnetic 
field with hardcore interaction [28], or the variant of the 
Calogero model considered in [29] . In any spatial dimen¬ 
sion, Richardson’s pairing model employing a contact 


pairing interaction potential allows for exact solutions 
l30l . Numerically, it has been shown using the harmonic 
interaction model 133(33 and its time-dependent gener¬ 
alization that the approach of (R-)MCTDHB can solve 
the generally time-dependent many-body problem to an 
in principle arbitrary degree of accuracy in ID (331 and 
also in 2D [ 31 ]. 

In what follows, we exploit the fact that conden¬ 
sate fragmentation is a genuine many-body phenomenon 
whose very existence relies on both first- and higher-order 
correlations. We show that the fact of the gas being 
no longer dilute, and hence Bogoliubov theory break¬ 
ing down, depends strongly on the power law of two- 
body interactions. In the present study, a system will 
therefore be considered dilute if it is accurately described 
by a mean-field theory because fragmentation is absent. 
On the other hand, a system is said to be not dilute 
(or strongly correlated), if fragmentation is present and 
mean-field theories are not applicable. We thus use the 
degree of fragmentation as an indicator of the “many- 
bodyness” of an interacting quantum gas: When frag¬ 
mentation becomes significant, the description of the sys¬ 
tem on an effective single-particle level, that is within 
mean-field theories such as the Bogoliubov theory, be¬ 
comes inapplicable. 

We will show that, indeed, the degree of fragmenta¬ 
tion is both small and essentially density-independent for 
contact interactions. However, the more long-ranged the 
interactions become, the more significantly the fragmen¬ 
tation degree depends on the density and exceeds the con¬ 
tact interaction value at sufficiently small densities. The 
degree of fragmentation is thus increasing when density 
of particles in the gas disk decreases. While this is a well- 
known fact in one dimension and for contact interactions, 
in two dimensions and for long-range interactions this re¬ 
sult is a nontrivial and novel manifestation of many-body 
physics. 

In addition, we make explicit that besides the density 
the total number of particles N is an independent pa¬ 
rameter by which to assess the mean-field behavior of 
the system. We demonstrate that in any dimension, and 
for any power law of interactions, the degree of fragmen¬ 
tation, and thus the degree to which a genuine quantum 
many-body nature of the interacting system is manifest, 
is reduced with increasing N. This markedly differs from 
the functional behavior of the degree of fragmentation 
upon varying p, which strongly depends on spatial di¬ 
mension and power law of interactions. 

II. MODEL AND METHODS 
A. Setup 

The iV-body Hamiltonian reads 

N N 

H = ^h(r i )+ Vint(ri,rj). (1) 

i— 1 i<j =1 
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Here, h{ri) = (— + Vext(fi)); is the one-body Hamil¬ 

tonian that contains the kinetic and potential energy of 
the particle with index i and Vj nt (r,;. r ? ) is the two-body 
interaction for a pair of particles at r t and rj. For rea¬ 
sons of computational convenience, we rescale the Hamil¬ 
tonian, by using a typical length scale of the system L, 

j-2 

by 2 mL 2 ■ This defines the units of length and energy we 
employ; all quantities in what follows are dimensionless. 

We consider a two-dimensional hard-wall disk potential 
specifically given by 

Kxt(r) = | 10Q0 r ~ a , ( 2 ) 


The areal density for N particles is thus given by p — -Ay, 
Note that most recent advances in engineering atom traps 
allow for hard-walled “box” traps, enabling the experi¬ 
mental study of the ground state of homogeneous Bose 
gases, in various spatial dimensions [551155] . 

In view of the convergence issues for a true S poten¬ 
tial in two dimensions mm, the contact interaction is 
modeled as a normalized Gaussian as follows 


Hint {R = |n - rj \) = g- 


exp - 


(-£) 


27TCT 2 


(3) 


with width a = 0.25 chosen, such that it reproduces the 
physics of a contact interaction [57]. The long-ranged 
interaction potential of two particles at ry and Vj , re¬ 
spectively, is modeled by the regularized expression 


Vint {R = | Vi -Vj |) 


9 

R n + A n ’ 


(4) 


where A(= 0.07) is the cutoff and n is the power law of 
the long-range interaction potential. 

Fragmentation is computed by calculating the natural 
orbital occupation numbers N t . which are the eigenvalues 
of the one-body density matrix 

P {1 \r'\r) = ^iV a 0V vo) ’*(r , )<)'i Aro) (r-), (5) 


with the eigenvectors (j)^ 0 ^ representing the natural or¬ 
bitals. The degree of fragmentation (for two modes) is 
then defined as [25] 


T = 1 - 


IJVr-JVal 

N 


( 6 ) 


where IV 1: 2 are the two largest natural orbital occupa¬ 
tion numbers, defined by Eq. |5]). We note that we have 
verified that in the parameter ranges under investigation 
below, the occupation of a possible third orbital remains 
negligibly small, so that performing a truncation after 
two modes in the field operator expansion is justified. 


B. Numerical Methods 

We use the multiconfigurational time-dependent 
Hartree method for Bosons (MCTDHB) [25] in its recent 


recursive implementation (R-MCTDHB) [IS] to compute 
the solution of the Schrodinger equation to an in princi¬ 
ple arbitrarily large degree of accuracy. The MCTDHB 
method uses the following ansatz for the wave function 


l*> = £Oi(t)|tf;i>. (7) 

n 


Here, IT) is expanded as a linear combination of 
time-dependent, fully variationally determined and sym¬ 
metrized many-body basis states, 


1 

irc;*) = -p= n mr ( g ) 

VIl"l»a!°=l 


where a vector notation was invoked for the occupa¬ 
tions, n = (ni,..., hm) T , which fulfill N = 

The many-body basis states |n; t) are obtained by apply¬ 
ing a symmetrization operator to Hartree products that 
are built from at most M distinct orthonormal orbitals 
{(j> a (r; t); a = 1,..., M}. With this ansatz one tackles the 
time-dependent many-body Schrodinger equation, using 
the time-dependent variational principle [39]. As the re¬ 
sult of the variation of the action functional, two sets of 
equations of motion, one non-linear integro-differential 
set for the time-evolution of the M single-particle ba¬ 
sis states or orbitals {<j> a (r\t);a = 1, ...,M} and one 
linear set for the time-evolution of the ( iV+ ^r /_1 ) coef¬ 
ficients {Cfj(t)}, are obtained, see Ref. [25] for a detailed 
discussion. The two sets are coupled, because the time- 
evolution of the orbitals depends on the matrix elements 
Pup(t), p a pis{t) which are functions of the coefficients, 
whereas the time-evolution of the coefficients depends 
on the matrix elements of the one-body and two-body 
Hamiltonian, h a p,W a pjs- The simultaneous and self- 
consistent solution of the system of MCTDHB equations 
of motion for both {Cfi(t)} and {<^ a (r;f)} is equivalent 
to solving the full (time-dependent) many-boson problem 
if convergence with respect to the number of variational 
parameters is achieved, see [Ml 133] • This is a conse¬ 
quence of the formal exactness of the ansatz in Eq. |t|; 
in the case of M —> oo, it covers the full iV-boson Hilbert 
space. Note that in the case of M = 1, Equation 0 
is a mean-field product state and therefore the set of 
MCTDHB equations boils down to the time-dependent 
Gross-Pitaevskii equation. 

The solution of the coupled set of equations for M > 
1 is a numerically demanding task for which the R- 
MCTDHB package [15] is used. 
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III. INTEGRABLE AND NONINTEGRABLE 
INTERACTIONS 

A. Two-dimensional disk geometry 

1. Fragmentation for various power laws of interaction 

We first remark, that for the present discussion, and 
especially for the presently considered mesoscopic sys¬ 
tems, N and p have distinct influence on the many-body 
state of the system. In the first-quantized form, N is a 
fundamental parameter of the many-body Hamiltonian 
while p is the ratio p = N/A\ here, the area A depends 
on the external potential. A change in N affects both the 
external potential and the interaction term while chang¬ 
ing A requires changing only the potential term and is 
independent of the interaction term. Since the N- and 
A-dependence of the system are different, the value of p 
does not uniquely determine the many-body wave func¬ 
tion, which also explicitly depends on the value of N. 

Hence, to analyze the density dependence of fragmen¬ 
tation independent from the influence of N, we keep 
N = 100 fixed and first vary the area A of the disk, 
by changing its radius a, in the 2D disk potential Eq. ([2]). 

The density dependence of fragmentation is assessed 
using a contact interaction and interactions with power 
laws n = 1,2,3 in Eq. Q. In principle, a comparison 
between different interaction potentials can be made with 
respect to the same coupling constant g. However, a more 
useful comparison can be obtained by using the following 
integrated measure of the 2D interaction potential V e g = 
2?r f V mt (R)RdR (oc g), since Vj nt enters the many-body 
Hamiltonian as the kernel of the integrals W a 0 7 s , see 
Ref. [22] and Eq. ([£]) below. The coupling constants are 
then renormalized such that they correspond to the same 
V e g for different interaction potentials. 

In Fig. [l] we compare the variation of fragmentation 
T with density p, holding V e g = 2n x 11.6 constant. 
For contact interaction (cf. Eq. ([3])), there is hardly a 
clearly discernible variation of J- with p. For long-ranged 
interactions (n = 1, 2 in Eq. Q), on the other hand, the 
degree of fragmentation decreases with p. For the 2D- 
integrable n = 3, and n > 3 (not shown), the influence 
of density on fragmentation is again minimal. 

This leads us to conclude that the influence of density 
on the degree of fragmentation for two spatial dimensions 
depends on the integrability of the interaction potential; 
for integrable interaction Eqs. and ([4]) for n > 3, the 
variation of T with p is minimal, while for nonintegrable 
interactions, n = 1,2 in Eq. Q, the fragmentation de¬ 
gree decreases with density. We also note that for n = 1, 
the most long-ranged case, the extent of the decrease 
of P(p) is significantly stronger than for n = 2, that is 
the more long ranged the interaction is, the stronger T 
varies with density. The n = 1 curve, while starting from 
higher values for small densities, crosses the curve for 
contact interactions, decreasing to a lower value of J- for 
higher densities. The n = 2 curve again starts with a 
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FIG. 1. (color online): The variation of fragmentation with 
2D density p = for N = 100 particles, comparing contact 
interaction and interactions of power law n = 1, 2, 3, at fixed 
V e f[ = 2n x 11.6. 
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higher value of J- for low density, decreasing slowly and 
approaches the contact interaction curve for large densi¬ 
ties. For n = 3, T(p) is highest for all densities. Finally, 
for sufficiently large densities, the n = 1 interaction is 
the least fragmented and the integrable n = 3 interac¬ 
tion is the most fragmented. We have verified that these 
qualitative features are similar for any number of parti¬ 
cles N and interaction coupling g with only the value of 
fragmentation J- being different. 


2. Effective Thomas-Fermi parameters 

The degree of fragmentation for a fixed N depends 
on the interplay between the kinetic energy and interac¬ 
tion energy. Thus, an important measure in determining 
the fragmentation of the system is the Thomas-Fermi pa¬ 
rameter represented by the interaction energy and kinetic 
energy, P TF = f^. 

Figure [2] shows the variation of the Thomas-Fermi pa¬ 
rameter Ptf with density p for fixed N = 100. One 
can observe a significant difference between integrable 
and nonintegrable interactions in particular in the low 
density regime. For contact interaction, there is a very 
slow decrease of Ptf with p throughout the entire den¬ 
sity regime; n = 3 shows similar behavior, with Ptf re¬ 
maining roughly constant with p. A strikingly different 
functional dependence is seen for n = 1. While the large 
p behavior is similar to that of contact interactions, for 
small p , there is a very sharp increase of Ptf with de¬ 
creasing density. A similar tendency (albeit to a smaller 
degree) is observed for n = 2. This strong increase of the 
Ptf for decreasing density, implying the stronger influ¬ 
ence of the interaction energy as compared to the kinetic 
energy, manifests itself in the increasing degree of frag- 
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FIG. 2. (color online): The variation of the Thomas-Fermi 
parameter -Ptf(p) with 2D density p for a fixed N = 100, 
for contact interaction and long-range interactions of degree 
n = 1,2,3, in two spatial dimensions. 


mentation with lowering density seen for nonintegrable 
interactions, cf. Fig.[T| 

To obtain a qualitative understanding of the above dis¬ 
cussed dependence on p, we note that in the mean-field 
regime for contact interactions, both the interaction and 
kinetic energy go linearly with density and Ptf does not 
vary with density. In a many-body context, the equally 
linear dependence of kinetic and interaction energy does 
not strictly hold and deviations are observed. In our 
case, the interaction energy goes slightly faster than the 
kinetic energy leading to the small gradual decrease of 
Ptf with density. The difference from the mean-field be¬ 
havior is, however, much more significant for n = 1 and 
n — 2, with the largest effect for n = 1. Here, especially 
for n = l, the variation of E- m t with p is not of a simple 
power law form and overall increases slower compared to 

in, leading to the rapid rise of Ptf at small densities. 
For example, when n = 1, E- mt ~ p 045 while E^ n ~ p 07 . 
Similarly, for n = 2, Ei nt ~ p 0 - 8 while E^ in ~ p °- 9 . 


3. Density profiles 



FIG. 3. (color online): 2D density distributions of the disk 
and the respective ID cut at y = 0 for Top: average density 
p = 7.95 (a = 2), for contact interactions, n = 1, n = 2, 
n = 3, from left to right. Bottom: average density p = 1.27 
(a = 5), at the same values of n. Darker and brighter color 
indicates lower and higher density, respectively; cuts at y=0 
are also shown. 


interaction for any radius of the disk. As a result, the 
equalization of the densities does not occur, and the rim 
of the disk continues to have higher densities, as clearly 
visible in Fig. [3] (bottom panel). 

This difference manifests itself in the occupation of the 
energy levels and thus the degree of fragmentation ob¬ 
served. The short range of integrable interactions leads 
to an essentially localized behavior and long-range cor¬ 
relations between the bosons are very weak. Noninte¬ 
grable interactions, on the other hand, imply longer- 
ranged correlations: The sphere of the influence of inter¬ 
actions involves the contribution of a much larger num¬ 
ber of bosons. A given boson thus interacts with all the 
bosons in the disk irrespective of the latter’s radius and 
hence feels a stronger effective interaction. In more for¬ 
mal terms, the kernel of the matrix elements of the two- 
body interaction, 


The difference between the influence of the integrable 
and nonintegrable interactions on fragmentation can be 
connected to the spatial density profile. For large densi¬ 
ties (small radius), the influence of the repulsive interac¬ 
tions for all values of n is felt on the entire disk, leading 
to a higher density in the rim of the disk (see Fig. [3] top 
panel). As the radius increases (thus lowering p ), the in¬ 
fluence of the short-ranged interactions (contact, n = 3) 
on the whole area of the disk is reduced, because the 
radius becomes much larger than any length scale asso¬ 
ciated to interactions, and the density hence is rendered 
homogeneous. The nonintegrable interactions (n = 1,2), 
on the other hand, retain the long-range influence of the 


W a p-,s= f d 2 rd 2 r'(f> a (r)(f>p(r')Vi n t(r — r')(f>^(r')(j)s(r) 

(9) 

has a larger support for smaller values of n (longer ranges 
of interaction). Therefore the integrals Wk sq i are effec¬ 
tively larger such that the system tends to minimize the 
energy by occupying several fij and hence fragments. 


B. Comparison with the one-dimensional case 

While in 2D, for contact interactions, we observe the 
absence of a change in fragmentation with density, in ID 
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FIG. 4. (color online): Variation of fragmentation with den¬ 
sity for N = 100 ID bosons with contact interactions and 
long-range interactions. Inset: Density dependence of frag¬ 
mentation on dimensionality at equal density for N = 100 ID 
and 2D bosons, interacting with unit strength g = 1. The ID 
density p 1D is scaled by \ to directly compare with the 2D 
range of densities p. 

there is a significant influence of the density, in line with 
what is known from exact diagonalization studies m- 
We compare a system of N = 100 bosons in the disk 
geometry for the 2D case and N = 100 boson in a ID 
“tube” with periodic boundary conditions. 

1. Long-range interactions 

Our simulations clearly reveal the difference between 
systems with long-range interactions in ID and 2D. In 
Fig. @ (main), we show the effect of the density p 1D = 
^ on fragmentation T in ID for a constant V e s = 
JV in t(R)dR = 5.15 with n. Fragmentation decreases 
with density for all n. For n = 1 in Eq. Q, the de¬ 
crease with density is most rapid. For n = 2, 3 (which 
are integrable in ID), the decrease is similar especially for 
larger densities. This is in complete accordance with our 
observations in the 2D bosonic gas: Concerning the de¬ 
pendence of T on p, integrable interactions demonstrate 
in either dimension similar behavior when compared to 
nonintegrable interactions. 

2. Contact interactions 

For contact interactions, in Fig. [4] (inset), we compare 
the variation of fragmentation with density for ID and 
2D gases for N = 100 and g = 1 (for the sake of a di¬ 
rect comparison with the 2D computations, we rescaled 
the ID density p 1D in the inset of Fig.Ulby a factor f). 
We observe a substantial decrease of .A for the ID case 



Number of particles N 


FIG. 5. (color online): (a) Number dependence of fragmenta¬ 
tion for contact and long-range interactions for 2D. (b) Num¬ 
ber dependence with fragmentation normalized to its value at 

N 0 = 10. 

with increasing density p 1D . The comparison to 2D with 
its absence of a significant variation of T with density 
again highlights the marginal nature of the two dimen¬ 
sional system. We note that in 3D, we anticipate that 
fragmentation is on one hand very small [40] , and on the 
other hand also decreases with decreasing density. 


IV. NUMBER DEPENDENCE OF 
FRAGMENTATION 

Fragmentation is a mesoscopic phenomenon in the 
sense that for fixed interaction coupling, increasing N 
always leads to a decrease of fragmentation. We have 
observed that this is a universal phenomenon, and true 
for both integrable and nonintegrable interactions, for the 
presently considered ground states of ID and 2D bosonic 
gases contained in a single trap. 

We illustrate in Fig.[5](a) how the fragmentation T 
varies in 2D with the number of particles N for contact 
and long-range interactions n = 1,2,3, cf. Eqs. 00' 
at a fixed V e s = 2n x 11.6. A decrease in fragmenta¬ 
tion is observed for all cases. Note that although with 
increasing N and constant coupling g , the mean-field in¬ 
teraction Vmf = (AT — l)p increases, this does not lead 
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to an increase of fragmentation T. Instead, we find that 
the particle number N dominates over the effects of all 
other parameters affecting the fragmentation T. That 
is, the variation of T with other parameters like inter¬ 
action strength g and long-rangedness n is subdominant 
in comparison. We stress that this qualitative behavior 
of the number dependence is universal in that it is true 
regardless of densities, interaction strength, and other 
parameters. Note that this fact is not captured by a 
(single-parameter) variational theory (4Tj , which demon¬ 
strates the sensitivity of fragmentation on the solutions 
of the many-body equations being fully self-consistent. 

In order to assess the relative decrease of fragmentation 
T, we plot in Fig.[5])b) the ^-dependence of T relative 
to the reference value set at N = 10. The contact inter¬ 
action curve shows the smallest decrease with N, closely 
followed by n = 1. The n = 2 and n = 3 show very sim¬ 
ilar dependence, and decay more rapidly with N when 
compared to n = 1. For constant interaction coupling, 
the decrease of the fragmentation degree is well approx¬ 
imated by T(N)/T(N 0 ) ex cf. Fig.lEFb). We finally 
mention that, when the effective mean-held interaction 
Vmf = (N — l)g were kept constant, then T falls off 
much faster (not shown), approximately as . 


V. CONCLUSION 

Using condensate fragmentation as a diagnostic tool, 
we have assessed the degree to which an interacting sys¬ 


tem of bosons can be considered to be dilute, and there¬ 
fore describable by mean-field theory and a single rnacro- 
scopically occupied orbital. We have confirmed that the 
essential independence of the diluteness parameter in a 
2D gas 0 is reflected by an identical behavior of the de¬ 
gree of fragmentation T in a bosonic quantum gas, for 
integrable interactions. However, for nonintegrable in¬ 
teractions, a significant dependence of T on the density, 
namely a significant decrease of the degree of fragmen¬ 
tation with increasing density obtains, and increasingly 
so for longer-ranged interactions. Due to the fact that 
fragmentation is a genuine many-body phenomenon, the 
degree of fragmentation thus represents a quantitative 
measure to which extent long-range interactions lead to 
strong correlations in two spatial dimensions. We have 
also found that in the large N limit, in a self-consistent 
approach, mean-field theory again becomes valid. Frag¬ 
mentation is therefore both a mesoscopic and intrinsically 
many-body phenomenon. Finally, the detection of the 
fragmentation we predict can be performed by measuring 
density-density correlations after time of flight expansion 
of the cloud [32] . 
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